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Abstract 

We report the result of a computation of the relation between the renor¬ 
malized coupling in the MS scheme of dimensional regularization and the bare 
coupling in the standard lattice formulation of the SU(Y) Yang-Mills theory 
to two-loop order of perturbation theory and discuss some of its implications. 
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1. The possibility to compute the running coupling in asymptotically free 
theories through numerical simulations of the corresponding lattice theories 
has recently attracted much attention and a lot of progress has been made 
in this held [1-16]. Different computational strategies are being pursued and 
several non-perturbative dehnitions of the coupling have been considered. In 
the pure SU(2) and SU(3) gauge theories the so far most precise and complete 
results have been obtained for a coupling Q!sf(?) related to the Schrodinger 
functional of the theory [2-5,8,9]. 

At sufficiently large momenta q any two running couplings can be ex¬ 
panded in powers of each other and the choices which one has made then 
become largely irrelevant. In particular, through the series 

"ms('^^) = «SF(g) + Ci(s) asF(g)^ + C2(s) asFiqf + • • • (1-1) 

(where s is an adjustable scale factor), one can convert from the SF scheme to 
the MS scheme of dimensional regularization [17,18], the currently most widely 
used scheme for the analysis of experimental data in high-energy physics. 

The one-loop coefficient ci(s) in the Yang-Mills theory with gauge group 
SU(2) and SU(3) has been calculated in refs.[2,5]. Considering the precision 
which can presently be reached in the numerical simulations, it is now desirable 
to extend the expansion (1.1) to two-loop order. In principle this could be 
done by evaluating the Schrodinger functional to this order using dimensional 
regularization. Inspection shows, however, that such a calculation would be 
technically delicate and computationally non-trivial. It seemed much safer to 
us to hrst expand the Schrodinger functional in the lattice theory in powers 
of the bare coupling oq and to combine the result with the series 

cif^(s/tt) = ciQ + di(s) ciq -|- d2(s) ciq , (1-2) 

which one derives on an inhnite lattice with spacing a. 

For the case of the SU(2) theory the expansion of the lattice Schrodinger 
functional to two-loop order has just been completed by Narayanan and Wolff 
[19,20]. For gauge group SU(Y) the one-loop coefficient di(s) appearing above 
has been computed nearly 15 years ago by A. and P. Hasenfratz [21] and an 
estimate for d 2 (s) has been obtained by Ellis and Martinelli [22-24] on the 
basis of “tadpole dominance” (as hrst proposed by Parisi [25]). We have now 
performed a complete calculation of d 2 (s) in the SU(Y) theory and discuss 
our result in this letter. The details of the computation will be published 
elsewhere. 
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2. We consider the standard SU(iV) lattice gauge theory on a four-dimen¬ 
sional hypercubic lattice with spacing a and Wilson action. The bare gauge 
coupling is denoted by go and oq = 6'o/4^- The Feynman rules in this theory 
can be deduced following well established lines (see e.g. ref.[26]). To determine 
the relation between oq and a set of correlations functions must be 

worked out to the desired order of perturbation theory and matched with the 
corresponding amplitudes which one obtains in the continuum theory using 
dimensional regularization. 

The analytical expressions for the ghost and gluon vertices are quite com¬ 
plicated on the lattice. As a result one tends to end up with Feynman inte¬ 
grands that easily hll several pages of paper even if all external Lorentz indices 
are contracted. Moreover, the Feynman parameter representation and other 
widely used tools to evaluate Feynman diagrams are of little help in the lattice 
theory, because the propagators are not simple quadratic functions in the loop 
and external momenta. It is, therefore, of crucial importance to arrange the 
computation in a most economical way. 

With this in mind we have decided to make use of the background held 
technique [27,28], which has previously been applied to compute the one- 
loop coefficient di(s) [29]. The idea of the method is to introduce a classical 
background gauge held and to couple it to the quantum helds in such a way 
that the system has extended symmetry properties. An important consequence 
of these symmetries is that the renormalization of the theory does not require 
any further counterterms besides those already required in the absence of the 
background held [28]. The background held should thus be regarded as an 
additional source held which allows one to probe the theory in a particularly 
interesting way. 

For our purposes the most suitable quantity to consider is the background 
held propagator, i.e. the second derivative of the effective action with respect 
to the background held at vanishing sources. To determine the matching of 
the gauge coupling in different schemes it is in fact sufficient to work out this 
amplitude to the desired order of perturbation theory. In addition the gauge 
held propagator must be calculated to lower orders to control the renormaliza¬ 
tion of the gauge hxing parameter. At two-loop order the computation then 
amounts to the evaluation of 31 propagator type Feynman diagrams on the 
lattice and a smaller number of diagrams in the continuum theory. Note that 
we do not set the external momentum to zero (in contrast to ref. [29]). The 
Feynman integrals are, therefore, perfectly well-dehned and no intermediate 
infrared regularization is required. 
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To evaluate the lattice Feynman integrals in the continuum limit, we 
have hrst extracted the dependence on the external momentum analytically 
and then computed the remaining parts of the integrals numerically through a 
newly developed position space technique. We expect this method to be more 
widely applicable and shall describe it in a separate publication. With little 
computer time the integrals which occurred could in this way be calculated to 
a numerical precision better than 7 signihcant decimal places. 

To guarantee the correctness of our results all lattice diagrams were cal¬ 
culated independently by each author, using different sets of algebraic and 
numerical programs on different computer systems. We checked algebraically 
that the quadratically divergent parts cancel in the sum of all diagrams and 
that the coefficients of the logarithmically divergent terms add up to the value 
expected from the Callan-Symanzik /3-function. The background held propa¬ 
gator in the MS scheme, which is also needed in our computations, has previ¬ 
ously been worked out to two-loop order by Ellis [24]. The formula quoted by 
Ellis completely agrees with our results. 

The outcome of our calculations is summarized by 

lliV TT 

=-ln(s)-TfciiV, (2.1) 

IV ; bTT ^ ^ 2iV ’ ^ ^ 

17/V^ 

d2(s) = di(sf - ^ 2^2 + h + ( 2 . 2 ) 

where 

h = 2.135730074078457(2), (2.3) 

k2 = -2.8626215972(6), (2.4) 

h = 1.24911585(3). (2.5) 

We do not comment on the size of d 2 (s) here, but would like to hrst put our 
result in the context of some concrete applications and shall then discuss the 
magnitude of the two-loop corrections in the relevant formulae. 

3. As already mentioned in sect. 1, we are mainly interested in the ex¬ 
pansion (1.1) relating the MS coupling to the SF coupling. By combining our 
results with those of Narayanan and Wolff [19,20], we can now calculate ci 
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and C 2 for N = 2 and obtain 


Ci(s) = -^ln(s) + 0.94327(4), 

OTT 

(3.1) 

C 2 (s) = ci(s)^ - ln(s) -7 0.5216(5). 

(3.2) 


For N = 3 the expansion of osF in powers of the bare coupling is currently 
only known to one-loop order so that in this section attention will be restricted 
to the SU(2) theory. 

Before making practical use of the series (f.f) the scale factor s should 
be hxed so that the low-order terms are reasonably small (as far as possible). 
If we choose 

s = Aj^/AsF = 2.24385(8), (3.3) 

one gets ci(s) = 0 and C 2 (s) = 0.0576(5), i.e. in this case the known pertur¬ 
bative corrections are indeed very small. Note incidentally that choosing s to 
be equal to the appropriate ratio of A-parameters proved to be successful in 
other instances as well, where the applicability of the perturbation expansion 
could be controlled through numerical simulations [8]. 

The SF coupling in the SU(2) theory has been accurately determined in 
ref.[8] over a large range of momenta given in units of Sommer’s scale tq [30]. 
Using the two-loop series (f.f) with s as given above, these results may now be 
converted to the MS scheme. At g = 200/ro for example (which corresponds 
to about 80 GeV in physical units), one obtains 

«Ms(?) = 0.1289(15)(3). (3.4) 

The hrst error here is the total statistical error as quoted in ref.[8], while 
the second is equal to Q!sf(?/'S)‘^ and thus indicates the expected order of 
magnitude of the three-loop correction in eq.(l.l). Since the two-loop term 
is anomalously small, it has practically no influence on the central value of 
a-^(q). The main impact of our computation of the two-loop correction is that 
the estimated error from the neglected higher order terms has been reduced 
from a few percent to a fraction of a percent. 

A remarkable property of the SF coupling is that its evolution in the 
momentum range covered by the numerical simulations is accurately described 
by the two-loop approximation to the Callan-Symanzik /3-function [3,8]. A 
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perfect fit of the data was in fact obtained by adding an effective three-loop 
term with coefficient 

= 0.35(12)/(47r)^ (3.5) 

(cf. appendix A for unexplained notation). From the above and the known 
value of 5^® [31] we may now calculate the exact coefficient 

5^^ = 0.1797(3)/(47r)^ (3.6) 

There is no reason to expect that 5®^ is equal to , but it is reassuring to see 
that both coefficients have the same sign and order of magnitude. Eq.(3.6) is 
certainly compatible with the numerical data. 

4. More direct ways of using the bare perturbation expansion (1.2) have 
recently been discussed in refs.[10,32]. The two-loop corrections in such ap¬ 
plications may now be calculated and we here consider two cases to illustrate 
this. For simplicity we set iV = 3 throughout this section. 

If we choose 

s = = 28.80934(1) (4.1) 

(as suggested by previous experience), the one-loop term in eq.(1.2) vanishes 
and d 2 {s) = 4.44915(1). One might think that with a different scale factor 
the rather large two-loop coefficient can perhaps be reduced signihcantly. But 
this is not the case, since d 2 {s) > 4.31303(1) for all s. 

In view of the size of the two-loop correction and the large scale factor 
(4.1) one may be hesitating to make practical use of this form of the expansion. 
The problem with the scale factor disappears if instead of oq one uses the 
“improved” bare coupling [25,32] 


do — cto/P (4.2) 

as the expansion parameter. To two-loop order the (normalized) plaquette 
expectation value P appearing in this dehnition is given by [33,34] 

^ dvr / , IGvr^A , , , 

P = 1 —4” ( ^^2 H g —) Oq 4“ • • • 5 (4-3) 

where k 2 is the constant (2.4) introduced earlier. By combining eq.(1.2) with 
eqs.(4.2),(4.3) the desired expansion of Q!j^(s/a) in powers of do is obtained 
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straightforwardly. We again choose the scale factor so that the one-loop term 
vanishes and end up with 

«ms(' 5/®) = «o + 2.18505(1) X s = 2.63285(1). (4.4) 

This expansion looks indeed more comfortable although the two-loop coeffi¬ 
cient is still rather large. In the relevant range of bare couplings, ~ 1, the 
correction is 3-4%. 

The introduction of the “improved” bare coupling (4.2) has partly been 
motivated by the observation that the dominant contributions to the one-loop 
coefficient di(s) come from tadpole diagrams and that these cancel when oq % 
replaced by do. Assuming tadpole dominance at two-loop order, as proposed 
by Ellis and Martinelli [22-24], however results in an estimate of the coefficient 
quoted in eq.(4.4) which is too small by more than a factor 3. An uncritical 
application of tadpole dominance can thus be rather misleading. 

5. The extension of our calculations to lattice QCD with Wilson quarks 
is straightforward. We also expect that the expansion of osF in powers of 
the bare lattice coupling can be worked out for this case to two-loop order 
along the lines explained in ref.[20]. Such computations would certainly be 
welcome, but it may take a while until the precision on osF which can be 
reached in numerical simulations of full QCD is comparable to the expected 
order of magnitude of the two-loop correction in eq.(l.l). 

We would like to thank Claus Vohwinkel for helpful discussions and Anton 
van de Ven for performing some additional checks on our computations. 

A. For any running coupling a = ^^/dvr depending on a momentum scale 
q we dehne the associated /3-function through 

/3(5) = ?|^- (A.l) 


In perturbation theory we have 


?4 = 0 

with 5o = llA/487r^ and hi = 17A^/3847r‘^ being the universal coefficients. All 
other (three-loop and higher order) coefficients depend on the chosen scheme. 
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The A-parameter is defined by 


A= lim (A.3) 

q^oo 

and if 

g’" =f + rig* + r2f + ... (A.4) 

is another running coupling, we have 

A'= Ae’'i/2fco^ (A.5) 

b'^ = b2 - biri+bo{r2 - rl). (A.6) 
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